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Abstract 

A remarkable feature of fluid dynamics is its relationship with classical dynamics and statistical 
mechanics. This has motivated in the past mathematical investigations concerning, in a special 
way, the "derivation" based on kinetic theory, and in particular the Boltzmann equation, of the 
incompressible Navier-Stokes equations (INSE). However, the connection determined in this way 
is usually merely asymptotic (i.e., it can be reached only for suitable limit functions) and therefore 
presents difficulties of its own. This feature has suggested the search of an alternative approach, 
based on the construction of a suitable inverse kinetic theory (IKT; Tessarotto et al., 2004-2008), 
which can avoid them. IKT, in fact, permits to achieve an exact representation of the fluid equations 
by identifying them with appropriate moment equations of a suitable (inverse) kinetic equation. 
The latter can be identified with a Liouville equation advancing in time a phase-space probability 
density function (PDF), in terms of which the complete set of fluid fields (prescribing the state of 
the fluid) are determined. In this paper we intend to investigate the mathematical properties of 
the underlying finite- dimensional phase-space classical dynamical system, denoted Navier-Stokes 
dynamical system, which can be established in this way. The result we intend to establish has 
fundamental implications both for the mathematical investigation of Navier-Stokes equations as 
well as for diverse consequences and applications in fluid dynamics and applied sciences. 

PACS numbers: 03.65.-w,05.20.Dd,05.20.-y 
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I. 1 - INTRODUCTION 



A fundamental theoretical issue in mathematical physics is the search of a possible finite- 
dimensional classical dynamical system - to be denoted as Navier-Stokes (NS) dynamical 
system - which uniquely advances in time the complete set of fluid fields which characterize 
a fluid system. In the case of an incompressible isothermal Newtonian fluid (also known as 
NS fluid) these are identified with the (mass) fluid velocity and the (non-negative) scalar 
fluid pressure, which in the domain of the fluid itself satisfy the incompressible NS equations 
(INSE). The reason why the determination of such a dynamical system is so important is 
that its existence is actually instrumental for the establishment of theorems of existence and 
uniqueness for the related initial-boundary value problem (INSE problem). In this paper 
we prove that, based on the inverse kinetic theory approach (IKT) developed by Tessarotto 
et al. (2004-2007 flEJEJ), the problem can actually be given a well-defined formulation. 
Main goal of the paper is the establishment of an equivalence theorem between the INSE 
problem and the NS dynamical system. Basic consequences and applications of this result 
are pointed out. In particular, contrary to the widespread view according to which the 
phase-space dynamical system characterizing the fluid fields should be infinite dimensional, 
here we intend to prove that a finite-dimensional classical dynamical system exists which 
advances in time the complete set of fluid fields. This is realized by the NS dynamical 
system. 



II. 2 - THE STRONG STOCHASTIC INSE PROBLEM 

For definiteness, let us assume that the complete set of fluid fields fluid fields {Z} = 
{p, V,p}, respectively denoting the mass density, the fluid velocity and the fluid pressure, 
describing an NS fluid, are local strong solutions of the equations 



P = Po, (1) 

V- V = 0, (2) 

AV = 0, (3) 

Z(r,t ) = Z (r), (4) 

z ( r J)\ d n = z ™( r ,t)\ d n' ( 5 ) 
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In particular, Eqs. (0Q)- (J5J) are respectively the incompressibility, isochoricity and Navier- 
Stokes equations and the initial and Dirichlet boundary conditions for {Z} , with {Z Q (r)} and 
{Z w (-,t)\ 9n } suitably prescribed initial and boundary-value fluid fields, defined respectively 
at the initial time t = t Q and on the boundary dfl. In the remainder, for definiteness, we 
shall require that: 

1. Q is coincides with the Euclidean space E 3 on R 3 (external domain) and dQ with the 
improper plane of R 3 ; 

2. / is identified with the real axis R (global domain) . 

We shall assume that the fluid fields are continuous in Q x I and fulfill the inequalities 

V > 0, (6) 

P > 0. (7) 
Here the notation as follows. iV is the NS nonlinear operator 

NV = ^V-F H , (8) 

with V and F# denoting respectively the Lagrangian fluid acceleration and the total force 
per unit mass 

^V=|V + V.W, (9) 
F H = Vp+— f + W 2 V, (10) 

Po Po 

while p and v > are the constant mass density and the constant kinematic viscosity. In 
particular, f is the volume force density acting on the fluid, namely which is assumed of the 
form 

f=-W(r,*) + fR, (11) 

4>(r,t) being a suitable scalar potential, so that the first two force terms [in Eq. ffTO]) ] can be 
represented as — Vp + f = — Vp r + f#, with 

p r (r,t)=p(r,t)-0(r,t), (12) 

denoting the reduced fluid pressure. As a consequence of Eqs. (pQ) and ([2]) it follows that the 
fluid pressure necessarily satisfies the Poisson equation 

V 2 p = S, (13) 



where the source term S reads 



S = - Po V ■ (V ■ VV) + V ■ f . (14) 

Here we shall assume, furthermore, that the fluid fields {Z} , together with the volume force 
density f and the initial and boundary conditions {Z Q (r)} , {Z w (r,t)\ dn } are all stochastic 
functions (see Appendix) of the form 

Z = Z(r,t,a), 
f= f (r, t, a) 

Z = Z (r,a) (15) 
Z w = Z w (r,i)\ d Q 

where a G V a are stochastic variables assumed independent of (r, t). Eqs.flTJ-flSJ) then denote 
the initial-boundary value problem for the stochastic incompressible Navier-Stokes equations 
{strong stochastic INSE problem). 



III. 3 - THE IKT-STATISTICAL MODEL 



A fundamental aspect of fluid dynamics lays in the construction of statistical models 
for the fluid equations 16) . In this connection a possible viewpoint is represented by the 
construction of the so-called IKT- statistical model able to yield as moments of the PDF the 
whole set of fluid fields {Z} which determine the fluid state and in which the same PDF 
satisfies a Liouville equation. Despite previous attempts (Vishik and Fursikov, 1988 [3( and 
Ruelle, 1989 ^|) the existence of such a dynamical system has remained for a long time an 
unsolved problem. 

This type of approach has actually been achieved for incompressible NS fluids [1]. Its 
applications and extensions are wide-ranging and concern in particular: incompressible ther- 

3019 



mofluids H, quantum hydrodynamic equations 
[10|], tracer-particle d yna mics for thermofluids 
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incompressible fluids 



12] , turbulence theory in Navier-Stokes fluids 
|l5| and applications of IKT to lattice-Boltzmann methods [lij . 



sco lb. '.-)\ ), phase-space Lagrangian dynamics 
17|, the evolution of the fluid pressure in 



and magnetofluids 



The IKT-statistical model is based on the introduction of a PDF depending on the local 
configuration-space vector r, /i(t) = /i(r, v,t, a) (1-point velocity PDF) defined on the 



phase-space I\ = Q x U [with U = M 3 ] and identified with a strictly positive function such 
that the complete set of fluid fields {Z} associate to the strong stochastic INSE problem 
can be represented in terms of the functionals (velocity moments) 

J dvG/i(t) = (l.V^aj.K^a)}, (16) 
u 

(Requirement #1; correspondence principle). Here G = {1, v,w 2 /2}while p\(r,t,a) > 
denotes the kinetic pressure 

pi(r, t, a) = p(r, t, a) + p (t, a) - 4>(r, t, a), (17) 

with po(t, a) > (the pseudo-pressure) a strictly positive, smooth, i.e., at least C^(I), real 
function and <f>(r,t, a) a suitably defined potential. In addition fi(t) is assumed to obey the 
Liouville equation - or inverse kinetic equation (IKE) according to the notation of Ref.j^j]) - 

L(v,v,t;f 1 (t))f 1 (t) = (18) 

(Requirement #2) with L(r, v, t; fi(t)) denoting the Liouville streaming operator 
L(r, v, t; fi(t))- = §- t ■ +^ ■ {X(x, t; fi(t))-} and F(fi(t)) a suitable smooth vector field 
defined in such a way that the moment equations of ( jTBl) obtained for G = {l,v,p w 2 /3} 
[u denoting the relative velocity u = u(r,t, a) = v — V(r, t,a)] coincide respectively with 
Eqs.()2]), ()3]) and again (|2]). This implies that the initial value problem associated to the 
vector field 

X(x,f;/i(f)) = {v,F(/i(t))}, (19) 



tlx 



X(x,t;A(t)) s 



x(t ) = x 

necessarily defines a dynamical system. In particular if x(t) = x(x ,t ,t) is the general 
solution of ( 12"U|) . this is identified with the flow 

T to ,t ■ x G -> x(t) = T to> iX = x(x , t , t) (21) 



generated by X(x, fi(t)). Furthermore it is assumed that Eq. (1181) admits as a particular 
solution fi(t) the Gaussian PDF 



where v^ h (r,t, a) = 2pi(r,t, a)/p . More precisely it is assumed that /m(x, t, a) is a particular 
solution of Eq.f JT8|) if and only if the fluid fields Z(r, t, a) are solutions of the strong stochastic 
INSE problem (Requirement #3). 

In the following we intend to investigate, in particular, the consequences of requirements 
(TTB|) . (ITS]) and (12"2"I) for the problem posed in this paper. 



IV. 4 - THE NAVIER-STOKES DYNAMICAL SYSTEM 

According to a certain misconception, dynamical systems for continuous fluids cannot be 
finite dimensional due to the fact that the fluid equations are PDEs for the relevant set of 
fluid fields {Z}. However, it is easy to show that this is not the case even in the so-called 
Lagrangian description of fluid dynamics. For a NS fluid this is realized by parametrizing 
the fluid velocity V in terms of the Lagrangian path (LP) r(t). In the present notation this 
is solution of the problem 

'f = V(r(t),t, a ), 
r(t ) = r , 

where is the Lagrangian derivative (J9]). As a consequence the NS equation becomes 

DVMt),t) „ , , s s 

^-Ai = F H (r(t),t), (24) 

which [with F#(r(t),t) considered prescribed] can be treated as an ODE and integrated 
along a LP yielding 

t 

V(r(t),t) = V(r ,t ) + J dt'F H (r(t'),t'). (25) 

to 

This permits to determine, for all r = r(t) and (r ,t ) £ fl x /, the vector field V(r,t). 
Therefore the finite- dimensional dynamical system (r ,£ ) — >■ (r(t),t) defined by Eq. (l23]) 
actually generates the time evolution of V(r, t, a). Nevertheless, in this description the fluid 
pressure is actually not directly determined [in fact this requires solving the Poisson equation 

We intend to show that a dynamical system advancing in time the complete set of fluid 
fields for a NS fluid is realized by the dynamical system T tojt [NS dynamical system; see 



Eq. (l2~T]) ]. Its precise definition depends manifestly on the vector field 
defining F(/i(i))] is achieved by the IKT approach developed in Ref. 



F(/i(t)). The task [of 
. As a consequence 
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it follows that F(/i(i)) can be non-uniquely determined 2|, |5[ in terms of a smooth vector 
field which is at least C^(Ti x I x V a ). 

In this fundamental mathematical result, we intend to prove here the equivalence 

between the strong stochastic INSE problem and the NS dynamical system. In other words 
the Liouville equation (lisp can be shown to be equivalent to 



fi(x,t,a) = J(t,a)f 1 (x(x,t,t ,a),t ,a), 
where x(t) = T <0)4 x is the general solution of (1201 for F = F(/i(t)) and where 



(26) 



J(t,a) = expij dt'— ^r— ■F(x(t',a),t',a;f 1 (t')) 



dx(t, a) 



<9x 



(27) 



, to 



is the Jacobian of the flow T to t [see Eq.( l2ip ]. Therefore, the NS dynamical system 
necessarily advances in time the PDF /m( x j t ,ct) so that it is identically a solution of 
the Liouville equation (ITS]) . The result can be established on general grounds, i.e., for an 
arbitrary vector field F((/i(i)) fulfilling Requirements #l-#3. The following result holds: 



THM.l - Equivalence theorem 

In validity of Requirements #1 — #3 the strong stochastic INSE problem is equivalent to 
the NS dynamical system WT\) . 
PROOF 

The proof is immediate. In fact, if 

x = x(t, a) = x(x o , t , t, a) (28) 

is the solution of Eq. (!20|) (which is assumed to exist and define at least a 
C^— diffeomorphism, its inverse transformation is simply 

x = x(x,t,t , a). (29) 

Therefore by differentiating Eq.f l26p it follows 

— /i(x,t,a) - — {J(t,a)fi(x(-x-,t,t ,a) 7 t ,a)} = 0, (30) 

which recovers Eq. (|26|) and admits as a particular solution /m(^) = /a/(x, t, a) when subject 
to the initial condition /m( x o, t Q ,a). Hence in terms of such an equation the NS dynamical 
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system advances in time the complete set of fluid fields. Therefore, the fluid velocity and 
the kinetic pressure at time t, i.e.,V(£) = V(r, t, a) and pi(t) = pi(r,t,a), follow from the 
moment equations (TTB"j) . Q.E.D. 



V. 5. CONCLUSIONS 

This work is motivated by the analogy between hydrodynamic description and the theory 
of classical dynamical systems. For greater generality the case of stochastic fluid equations 
has been considered. The problem of the equivalence between the initial-boundary value 
problem for incompressible Navier-Stokes equations and the Navier-Stokes dynamical system 
introduced in Ref.[2j has been investigated. Indeed, the theory here developed applies both 
to deterministic and stochastic fluid fields. In fact, in both cases the time evolution of fi 
is determined by a Liouville equation [see Eq. (fr8l) ] which evolves in time also the complete 
set of fluid fields (all represented in terms of moments of the same PDF). Contrary to the 
misconception according to which the phase-space dynamical system characterizing the fluid 
fields {Z} of a continuous fluid system should be infinite dimensional, here we have proven 
that the finite-dimensional NS classical dynamical system advances in time the complete set 
of fluid fields, determined in terms of velocity moments of the 1-point PDF /i(x, t,a). The 
;heory here developed applies generally to stochastic fluid equations. As shown elsewhere 
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16| . this represents a convenient treatment for the statistical theory of turbulence, 



historically referred to the work of Kolmogorov (Kolmogorov, 1941 [18|) and Hopf (Hopf, 
1950/51 [lj). 

The theory has important consequences which concern fundamental aspects of fluid dy- 
namics: 

• determination of the NS dynamical system advancing in time the complete set of fluid 
fields of a turbulent NS fluid [16]; 



construction of the initial conditions for the 1-point PDF fi 20]; 



determination of the time-evolution of passive scalar and tensor fields 



12|; 



construction of the exact equations of motion for ideal tracer-particle dynamics in a 
turbulent NS fluids [17| : 
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construction of multi-point PDFs for turbulent NS fluids [21] ; 



statistical treatment of homogeneous, isotropic and stationary turbulence based on 
IKT [22J. 



VI. APPENDIX - STOCHASTIC VARIABLES 

Let (S, E, P) be a probability space; a measurable function a :S — > V a , where V a C R k , 
is called stochastic (or random) variable. 

A stochastic variable a is called continuous if it is endowed with a stochastic model 
{g a ,Va} j namely a real function g a (called as stochastic PDF) defined on the set V a and 
such that: 

1) g a is measurable, non-negative, and of the form 

9a = 9a(r,t, ■); (31) 

2) if A C V Q is an arbitrary Borelian subset of V a (written A G B(V a )), the integral 

P a {A) = J dxg a (v,t,x) (32) 

A 

exists and is the probability that a Gv4; in particular, since a E V a , g a admits the normal- 
ization 

J dxg a (r, t, x) = P a {V a ) = 1. (33) 

v a 

The set function P a : B(V a ) — >■ [0, 1] defined by (132]) is a probability measure and is called 
distribution (or law) of a. Consequently, if a function f:V a — > Vf C ]R m is measurable, / 
is a stochastic variable too. 

Finally define the stochastic-averaging operator (-} a (see also [l^]) as 

(/)« = (/(y,0> a = f rfx^(r,t,x)/(y,x), (34) 
for any P^-integrable function /(y, •) : V^, — >■ M, where the vector y is some parameter. 
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